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O ' ABSTRACT 

A Chandra spectrum of the moderately active nearby F6 V star ir 3 Ori is used 
to study the coronal properties of mid-F dwarfs. We find that 7r 3 Ori's coronal 
emission measure distribution is very similar to those of moderately active G 
and K dwarfs, with an emission measure peak near logT = 6.6 seeming to be 
ubiquitous for such stars. In contrast to coronal temperature, coronal abundances 
C/3 \ are known to depend on spectral type for main sequence stars. Based on this 

r£3 [ previously known relation, we expected 7r 3 Ori's corona to exhibit an extremely 

strong "FIP effect," a phenomenon first identified on the Sun where elements 
2 | with low "First Ionization Potential" (FIP) are enhanced in the corona. We 

instead find that 7r 3 Ori's corona exhibits a FIP effect essentially identical to 
that of the Sun and other early G dwarfs, perhaps indicating that the increase in 
FIP bias towards earlier spectral types stops or at least slows for F stars. We find 
that 7r 3 Ori's coronal characteristics are significantly different from two previously 
studied mid-F stars, Procyon (F5 IV-V) and r Boo (F7 V). We believe 7r 3 Ori is 
more representative of the coronal characteristics of mid-F dwarfs, with Procyon 
being different because of luminosity class, and r Boo being different because of 
the effects of one of two close companions, one stellar (r Boo B: M2 V) and one 
planetary. 



Subject headings: stars: individual (it 3 Ori) — stars: coronae — stars: late-type 



— X-rays: stars 



INTRODUCTION 



Although details about stellar magnetic field generation remain uncertain, it is accepted 
that coronal magnetic fields are generated by a magnetic dynamo operating within the 
convection zone. For main sequence stars, the convection zone narrows towards earlier 
spectral types, di sappearing entirely a round a spectral type of A5 V. At this point, coronae 
disappear as well ( jSchmitt et al.lll985f ). emphasizing the central role that the convection zone 
dynamo plays in generating stellar coronae. 
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One general empirical goal of stellar coronal X-ray observations is to study how coronae 
change character as the convection zone narrows. Any clear correlation between coronal 
properties and spectral type tells us something useful about how the dynamo and the fields 
it generates change with spectral type, thereby providing potentially crucial information for 
theoretical models of how stellar magnetic dynamos operate in general. 

Coronal abundances represent one coronal property that appears to be tightly correlated 
with spectral type, at least for main sequence stars. On the Sun, it has long been known 
that coronal abundances are different from photospheric abundances, elements with low First 
Ionization Potential (FIP) having enhanced abund ances in the solar corona and solar wind 
( Ivon Steiger et al.l Il995t iFeldman fc Laming! 120001 ) . Analyses of X-ray and EUV spectra 



have found a similar "F I P effect" for some stars of low to moderate a ctivity ( ILaming et al. 



19961 ; iDrake et al.l 119971 ; ILaming &: Drake! Il999l ; iTelleschi et al.l 120051 ) . However, for many 



other stars, particularly extremely active ones, the FIP effect is either absent, or sometimes 
an inverse FIP effect is observed, where low-FIP elements have coronal abundances that 
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Initially, it appeared that coronal abundances were primarily dependent on activity, as 
the most active stars clearly tend to have inverse FIP effects instead of a solar-like FIP 
effect. T his activity de pendence of abundances for active stars is illustrated and discussed in 
detail by iGiidell ( 120041 . see Figure 37). The abundance variations are not simply a uniform 
variation of low-FIP elements relative to high-FIP ones. The high-FIP elements appear to 
fractionate with respect to each other as well, as illustrated by measurements of Ne/O. The 
yast majority of spect roscopically observed stellar coronae exhibit high ratios of Ne/O ~ 0.4 
( IDrake fc Testa! 120051). but relatively in active stars (e.g., the Sun, a Cen AB, and Procyon) 
show Ne/O « 0.2 jRobrade et al.l l2008h . 



As more abundance measurements have bee n made, it has become apparent that abun 



danc es are also highly dependent on spectral type (jGudelll2007l ; lGudel et al. 



2007; Telleschi et al 



20071 ). In order to explore this further we have experimented with a high-FIP/low-FIP abun- 
dance ratio, F}, ias , which quantifies a corona's FIP bias as the average abundance of four 
high-FIP elements (C, N, O, and Ne) relative to Fe; normalized to photospheric abundance 
ratios. For normal main sequence stars with coronal X-ray luminosities of logLx < 29.1 (in 
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ergs s " 1 ), there is a sur prisingly tight spectral type dependence for Fb ias ( IWood &: Linsky 
20101 ; IWood et al.ll2012l ). In this correlation, early G stars all have a solar-like FIP effect, 
which decreases towards late G and early K stars, reaching no FIP effect at all at a spectral 
type of K5 V. Later than K5 V, inverse FIP effects are observed, with the magnitude of the 
effect increasing into the M dwarfs. We will refer to this as the "FIP-Bias/Spectral-Type" 
(FBST) relation. At least for stars with logLx < 29.1, any activity dependence for Ff, ias is 
lost in the scatter of the measurements. 

The FBST relation would predict that the solar-like FIP effect should become even 
stronger as you move from solar-like early G stars to F stars. More pr ecisely, instead of a 
facto r of 4 enhancement of low FIP elements in the corona as for the Sun (IFeldman fc Lamingl 
2000), the FBST relation would predict a factor of 6 enhancement for mid-F stars (IWood &: Linsky 
20101 ). However, observations of two mid-F dwarfs seem to contradict this expectation. The 
first of these is Procyon (F5 IV- V), the nearest and most frequently observed F star, which 
has been observed spectroscopically by EUVE, Chandra, and XMM. All of these spectra 
demonstrate that Procyon's corona exhibits no significant FIP bias at all. It s coronal abun 



danc es are consistent with its photosp heric abundances (jDrake et al.l Il995t iRaassen et al. 



20021 ). More recently, using XMM data lMaggio et al.l (120111 ) have also found little difference 
in the abundances of low-FIP and high-FIP elements in the corona of r Boo (F7 V), although 
they suggest that both are depleted in the corona by similar amounts. 

The Procyon and r Boo measurements would seem to suggest that the FBST relation 
falls apart for the F dwarfs, possibly an effect of the thinning convection zones of such stars. 
This would be a very interesting result, but there are potential issues with both Procyon 
and r Boo regarding whether either can be considered representative of mid-F dwarfs. For 
Procyon, the problem is simply that it is not a pure main sequence star, as its "IV- V" 
luminosity class sug gests. Its radius of 2 .05 R is certainly larger than that of a true mid-F 
main sequence star (IKervella et al.l 12004 ) . 



For r Boo, the problem is the presence of two very different companions, one stellar and 
one planetary. The stellar companion, r Boo B, is an M2 V star that is unresolved from the 
primary in the XMM observation, which could in principle be contributing X- ray photons 
to th e spectrum. Given that M dwarfs are expected to have inverse FIP effects (IWood et al. 
20121 ) . contamination from r Boo B would be expected to weaken any FIP effect in the 
combined X-ray spectrum. As for the planetary companion, it is a very well-studied close-in 
giant planet, r Boo b, with an orbital period of only 3.31 days. There is alread y evidence that 
the presence of this mass ive planet is producing coronal activity on the star (jShkolnik et al. 
20081 ; IWalker et al.l 120081 ) . Is it possible that the planet could somehow be suppressing the 
expected strong coronal FIP effect as well? This is an exciting possibility, which makes it 
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all the more important to know what the corona of a truly normal mid-F dwarf looks like. 



To that end, we here present an analysis of recent Chandra observations o f the nearby 



F6 V star, tt 3 Ori, which lies at a distance of 8.1 pc f lPerryman et al.l Il997l ). The star 



has a rapid rotation rate of y sm i = 17.3 km s -1 , corresponding to P rot / sinz = 4.1 days 
( lAmmler-von Eiff fc Reinersll2012l ). Not surprisingly, given the known connection between 
stellar rotation and activity, t t 3 Ori also has a relatively high X-ray luminosity of log Lx = 
28.96 ( jSchmitt fc Liefkd l2004j ) . just barely within the logLx < 29.1 activity regime men- 
tioned above where the FBST relation remains tight. This is in fact one of the ten most 
luminous coronal X-ray sources within 10 pc, so it is curious that tt 3 Ori has never before 
been observed spectroscopically in X-rays. The tc 3 Ori spectrum promises to resolve the issue 
of what the coronal characteristics of moderately active mid-F main sequence stars really 
are. 



2. DATA ANALYSIS 

2.1. Observations of ir 3 Ori 

The Chandra X-ray telescope observed 7r 3 Ori twice in late 2010, the first time on 
2010 November 11 for 57.8 ksec, and a second time on 2010 November 23 for 20.1 ksec. 
The observations were performed with Chandra's LETGS configuration, combining the Low 
Energy Transmission Grating (LETG) and the component of the High Resolution Camera 
detector used for spectroscopy (HRC-S). 

The data were processed using version 4.3 of the Chandra team's CIAO software. The 
initial HRC-S image produced by the data processing consists of a zeroth-order image of the 
star, and an X-ray spectrum dispersed on both sides of the image, i.e., a plus and minus 
order spectrum. Further processing extracts one-dimensional spectra from the image, and 
then the plus and minus order spectra are coadded into a final spectrum covering the full 
5 — 175 A wavelength region of LETGS. 

Photons in the time-tagged zeroth-order image are useful for assessing source variability. 
We find no significant variability at all for tt 3 Ori, and the X-ray fluxes appear to be identical 
at the two observing times, so we simply coadd the spectra taken at the two times into a 
final spectrum representing a combined exposure time of 77.9 ksec. This spectrum is shown 
in Figure 1. 
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2.2. Line Identification and Measurement 



We u se version 6.1 of t he CH IANTI atomic database to identify emission lines in our 
spectrum (IDere et al.l 119971 . 120091 ). The identified lines are listed in Table 1, and identifica- 
tions are also provided in Figure 1. Counts for the detected lines are measured by direct 
integration from the spectrum, and these raw counts are then converted to photon fluxes. 
Both the counts and fluxes of the lines are provided in Table 1, together with la uncertain- 
ties. The third column of the table lists line formati o n tem peratures, based on ionization 
equilibrium computations from lArnaud fc Rothenffugi (119851 ) . 



There are many lines in Table 1 with large quoted uncertainties that represent ques- 
tionable detections. The emission measure (EM) analysis described in section 2.3 is used to 
confirm a posteriori that these marginal detections are at least plausible. There are many 
line blends in Table 1. The EM analysis and the line strengths in the CHIANTI database 
are used to assess what lines are contributing significant flux to each emission feature. The 
line blends involving different species are ultimately ignored in the final EM analysis. 



There are many density sensitive He-like triplets in the ir 3 Ori spectrum (INess et al. 



20021 ) . but only the O VII triplet near 21.8 A is sufficiently free of blends and with high 
enough signal-to-noise (S/N) for a decent density measurement. Unfortunately, even for 
O VII we find we can only quo te an upper limit for the electron density. Based on the 
collisional equilibrium models of iPorquet et al.l ( 1200 ll ). the A22.101/A21.807 line flux ratio 
corresponds to an electron density (in cm -3 ) of logn e < 10.70. Similarly low coronal densities 
have been found for many moderately active G and K dwarfs in the past ( IWood fc Linsky 
20061 . 120101 ). so there is no evidence that our moderately active F dwarf is any different from 
later type stars in terms of coronal density. 



2.3. Emission Measure Analysis 



Inferring coronal abundances and temperature distributions from line flux measurements 
requires an EM analysis, which is a nontrivial inve rsion problem. For this c alcula tion, we use 
version 2.6 of the PINTofALE software package ( iKashyap fc Drake! Il998l 120001 ). which was 
developed to solve the EM inversion problem using a Mar kov chain Monte Carlo approach. 
The analysis also makes use of CHIANTI line emissivities (IDere et al.lll997l 120091). Our EM 
analys is closely mirrors previous work, as described in particular detail by IWood fc Linsky 
( 120061 ). so we refer the reader to that paper for details about this computation. 



The analysis considers the effects of interstellar absorption on the observed line fluxes. 
There is no measured H I column density towards 7r 3 Ori itself, but there is a measurement 
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of log Nh = 17.93 (in cm -2 ) towards \ l Ori (IWood et al.ll2005l ). which is roughly at the 
same distance and in the same direction as 7r 3 Ori, so this is the value we assume. We note, 
however, that this value is low enough that interstellar absorption does not have a large 
effect on line fluxes. 

The line fluxes by themselves allow the computation of the shape of the emission mea- 
sure distribution, and the computation of relative coronal abundances. For the relative 
abundances, the natural reference element to use is Fe, considering the large number of Fe 
lines available in the LETGS spectrum. In the initial EM analysis we simpl y assume a so 



lar p hotospheric abundance for our reference element, log [Fe/H] = —4.50 (lAsplund et al. 



20091 ). Table 2 lists the coronal abundances relative to Fe computed in the EM analysis, 
and Figure 2 shows the derived emission measure distribution. The 90% confidence intervals 
shown for the emission measures and abundances in Figure 2 and Table 2, respectively, are 
indicative of the uncertainties in solving the inversion problem, based on the magnitude of 
the uncertainties in the measured line fluxes. 

Properly normalizing the EM distribution requires the measurement of an absolute 
coronal [Fe/H] ratio, which also allows the relative abundances in Table 2 to be converted to 
absolute abundances. Measurement of [Fe/H] requires an assessment of the line-to-continuum 
ratio in the LETGS spectrum, which can be done after the initial EM analysis is completed. 
Emission measures scale inversely with the assumed [Fe/H], since higher metal abundances 
mean less emission measure is required to account for the line fluxes. Thus, increasing [Fe/H] 
by a factor of two from log [Fe/H] = —4.50 to log [Fe/H] = —4.20, for example, uniformly 
lowers the EM distribution by a factor of two. 

Synthetic line-plus-continuum spectra can be computed from the EM distribution for 
various values of [Fe/H] and compared with the data, in order to see which leads to the 
best match to the observed continuum level. This is illustrated in Figure 3, which shows 
the best fit for log [Fe/H] = —4.32. This synthetic spectrum is also compared with the 
data in Figure 1. In LETGS data, higher order spectra are superimposed onto the first order 
spectrum, so the synthetic spectrum accounts for the higher order spectra as well, specifically 
orders 2-5. The higher order contributions are explicitly plotted in both Figures 1 and 3. The 
EM distribution in Figure 2 is the final normalized distribution assuming log [Fe/H] = —4.32. 



2.4. The Shape of the Emission Measure Distribution 



The most striking feature of the tc 3 Ori EM distribution in Figure 2 is a sharp peak 
at logT = 6.6. This peak actually seems to be very common for moderately active stars 
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like 7r 3 Ori. Many examples are provided by I Wood &: Linskyl (120061 |2010| ). but the best of 
these are £ Boo A (G8 V) and e Eri (K2 V). The EM distributions for these 2 stars are also 
shown in Figure 2. The £ Boo A EM distribution is a particularly good match for that of 
7r 3 Ori. This is perhaps not surprising considering that £ Boo A and 7r 3 Ori have similar 
X-ray luminosities of logLx = 28.86 and logLx = 28.96, respectively (ISchmitt fc Liefke 
20041 ). although the spectral types of the two stars are quite different. Thus, there is no 
evidence that the thinner convection zone of 7r 3 Ori has led to coronal temperatures any 
different from those of similarly active later type stars. 

Int erestingly enough, the log T = 6.6 peaks are also commonly observed for solar active 
regions (IWarren et al.ll201ll . |2012| ). In Figure 2, the three stellar EM distributions are com- 
pare d with a solar one, w hich represents the average of 15 active region EM distributions 
from I Warren et al.l (120121 ). who emphasized the prevalence of logT = 6.6 EM peaks on the 
Sun. The shape of the solar active region EM agrees beautifully with those of 7r 3 Ori and 
£ Boo A for logT < 6.7, a lthough the three s tellar EM distributions in Figure 2 have much 
higher EM at logT > 6.7. iDrake et al.l (120001 ) have previously noted the similarity between 
the EM distributions of solar active regions and those of £ Boo A and e Eri, and have sug- 
gested that it implies that moderately active stars like these are covered almost completely 
with solar-like active regions. 

Is it possible that these seemingly ubiquitous logT = 6.6 peaks are somehow a system- 
atic artifact of the PINTofALE computations that produced them? Not only was PINTo- 
fALE's Markov chain Monte Ca rlo approach used to compu te the stellar EM distributions in 
Figure 2, but it was also used by I Warren et al.l (120111 . 120121 ) to compute the solar ones. Note, 
however, that I Warren et al.l (120111 . 120121 ) analyzed a completely different spectral region and 
a different set of emission lines than those used in the stellar analyses, suggesting that the 
logT = 6.6 peaks are not just a product of working with a particular set of lines. We explore 
this issue further using simulations. These simulations involve taking a known EM distribu- 
tion, computing line fluxes from that distribution for the set of lines in Table 2, and then 
plugging those line fluxes back into PINTofALE to see if the EM inversion procedure can 
recover the original EM distribution. Uncertainties in the line fluxes are assumed to be the 
same as the relative uncertainties measured for the it 3 Ori lines in Table 1. Some results of 
these simulations are shown in Figure 4. 

The dotted line in Figure 4's "Simulation 1" is simply the 7r 3 Ori distribution from 
Figure 2. The solid line shows the result when that is used as the input EM distribution for 
the simulation. The computed EM distribution agrees encouragingly well with the original 
distribution, being within the displayed error bars for 15 of the 17 temperature bins. Most 
importantly, the simulation demonstrates PINTofALE's ability to recover the sharp log T = 
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6.6 EM peak. In Simulations 2 and 3 we simply shift the peak to logT = 6.5 and logT = 6.7, 
respectively, and in Simulation 4 we remove the peak entirely, assuming equal EM values in 
the logT = 6.5 — 6.7 range. Although the reconstructed EM values are still generally within 
error bars, we do note that for Simulations 2-4 the EM is systematically overpredicted at 
logT = 6.6 and underpredicted at logT = 6.5 and logT = 6.7. 

We conclude, therefore, that there may indeed be a modest tendency for PINTofALE 
to overemphasize the height and narrowness of the logT = 6.6 peaks, though the magnitude 
of this effect is not necessarily inconsistent with the inferred 90% confidence intervals. We 
can only speculate as to the source of this tendency. It does so happen that logT = 6.6 
is the central formation temperature of the strongest lines in the LETGS spectrum, the 
Fe XVII lines between 15-17 A, so perhaps the need to fit those line fluxes particularly 
well leads to the modest overemphasis of EM at logT = 6.6, to the detriment of EM at 
logT = 6.5 and logT = 6.7. In any case, despite evidence for this small systematic error, the 
simulations provide support for the existence of a distinct maximum of EM near logT = 6.6 
for moderately active stars and solar active regions. 

The question then becomes, why is logT = 6.6 such a popular temperature for both a 
variety of solar active regions, and for moderately active stars of various spectral types? Is 
there something special about this particular temperature? One notable thing a bout it is that 
there is a local minimum in the radiative cooling curve near logT = 6.6 (e.g., IColgan et al. 
20081 ; iDudfk et al.ll201ll ). The cooling curve relates how emissivity varies with temperature. 
Plasma is stable against radiatively driven thermal instabilities where the cooling curve has a 
shallower slope than the temperature dependence of the heating function. The cooling curve 
has a steep negative slope for a variety of abundance sets in the range logT = 6.1 — 6.5, 
which could imply thermal instability unless the heating curve has an even steeper negative 
slope at those temperatures. Thus, logT = 6.6 is possibly the coolest stable temperature 
above logT = 6.1. Plasma initially between these two temperatures will most likely either 
heat up or cool down until a stable temperature is reached. 

Perhaps the ubiquity of the logT = 6.6 peak implies something about the commonality 
of active region loops on the Sun and solar-like main sequence stars, consistent with the idea 
that coronal emission from moderately active stars like 7r 3 Ori and £ Boo A is d ominated 
by solar-like active regions that largely cover the stellar surface ( Drake et al. 200oh . On the 



Sun, such loops are often found to have densities of n e ~ 10 10 cm 3 (e.g., IWinebarger et al. 
201ll ). Such densities are consistent with those commonly measured in stellar coronae as well, 
at least at temperatures of logT < 6.7 JNess et al.lbooi IWood k Linskvl bood . boioh . An 
explanation in terms of thermal stability presupposes that coronal plasma exists in structures 
(loops) that can support a quasi-steady state. No such explanation would be possible in 
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a corona heated and supplied by more dynamic phenomena such as Type II spicules (e.g . 



de Pontieu et al.ll201ll ; Mcintosh et allbOlll Martinez- Svkora et al.ll201ll ; Ijudee et al.ll2012l ) . 



In simple static loop mod els there is a relation between a coronal loop's length and its ape x 
temperature and density (IRosner et al.lll978l ; IVesecky et al.lll979l ; |Jordanlll980l ; lRealdl2010l ). 
Thus, if coronal loops in various solar and stellar contexts tend to have similar densities and 
lengths, then presumably they would have similar temperatures as well. This presumes, of 



course, that these s imple static models apply, which is questionable (e.g., iLenz et al.l Il999 



Schmelz et all 120091 ) . 



2.5. The Coronal Abundances of n 3 Ori 



As discussed in section 2.3, the EM analysis provides coronal abundances measured rel- 
ative to Fe for many elements. All of these measurements, in the form log [X/Fe], are listed 
in the second column of Table 2. We are interested in comparing coronal abundances with 
those of the photosphere. Thus, in the third column of Table 2, we list stellar photospheric 
abundance ratios relative t o Fe, in the form log [X/Fe] The photospheric abundances of 
Allende Prieto et al.l ( 120041 ) are measured using line-by-line comparisons of solar and stellar 
spectral features, so the abundance measurements are intrinsically relative to solar abun- 
dances. Thus, computing the stellar photospheric abundance ratios i n Table 2 requi r es the 
assumption of solar abundances. The assumed solar abundances of lAsplund et al.l (120091 ) 
are listed in the last column of Table 2. For N and S, there are no stellar photospheric 
measurements available. We assume these elements behave like high-FIP elements, such as 
O. Thus, we assume log [N*/N ] = log [S„/S ] = log [O*/O ] = 0.16. 

Neon is a special case. The fundamental problem with Ne is that even for the Sun there 
is no real photospheric abundance measurement, because there are no solar photospheric 
Ne absorption lines. Reference solar abundances for Ne, such as the one listed in Table 2, 
are in reality based on measurements from the solar corona and transition region (e.g., 
Schmelz et al.ll2005l ). We know very well that abundances can be di fferent f r om ph otospheric 
in the corona, but not necessarily in the transition region, where lYoungi (l2005aU bl) finds a 
low Ne/O abundance, and no evidence for a FIP effect in other elements. Nevertheless, 
we do know that high-FIP elements can be frac tionated in the solar at mosphere, as He is 
depleted by a factor of 2 in the slow solar wind (Ivon Steiger et al.lll995l ). This depletion of 
He (and also possibly of Ne) is also a feature of models for the solar FIP effect based on the 
ponderomotive force (jLamingj|2009l . 120121 ; iRakowski &: Lamina |2012| ). 



Complicating the issue further is the knowledge that moderately active and very active 
stars seem to univerally have coronal Ne/O ratios much higher than that of the Sun and other 
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similarly inactive stars ( IDrake fe Testa! 120051 ; IWood fe Linskyl 120061 ; iRobrade et al.l 120081 ) . 
raising the possibility that these stellar coronal Ne abundances may be a better measure 
of the true cosmic Ne abundance than the solar coronal Ne abundance usually referenced 
as such. We here follow the same convention as in ou r past papers a n d ass ume that the 
average stellar coronal Ne/O = 0.41 abundance ratio of IDrake fe Testal (120051 ) is indicative 
of the true cosmic abundance of Ne, and therefore applies to solar and stellar photospheres 
as well, including for 7r 3 Ori. Thus, the logpNe/Fe]^ ratio in Table 2 is derived assuming 
Ne/O = 0.41, instead of the much lower solar ratio of Ne/O = 0.17. The validity of this 
assumption is very much debatable, but it does lead to coronal Ne abundances in better 
agreement with those of other high-FIP elements (see Figure 5). 

After subtracting the logarithmic photospheric abundance ratios, the logarithmic coro- 
nal abundance ratios in Table 2 are in Figure 5a plotted versus FIP. The low-FIP elements 
are all roughly consistent with the reference low-FIP element, Fe. The high-FIP element 
abundances are lower by about —0.6 dex, consistent with a FIP effect nearly identical to 
that of the Sun. The average value of the four high-FIP elements in Figure 5a is plotted 
explicitly as a horizontal dot-dashed line in the figure. This is the "FIP bias" quantity, F bias , 
that we have used in the past to reduce coronal abundanc es to a single number, which is then 
used to study the FBST relation described in section 1 ( IWood fe Linskyl l2010t IWood et al. 
2012h . For 7T 3 Ori, F bias = -0.55. 



The F}, ias quantity has no clear variation wi th activity for main sequence stars with 
logL x < 29.1 (see Figure 4 in IWood et al.l 120121 ). but it is worth noting that of the four 
high-FIP elements that go into computing Fu as (C, N, O, and Ne), Ne is once again an ex- 
ception, as we know that the coronal Ne/O ratio does vary within the logLx < 29.1 sample, 
suggesting activity-dependent fractionation of Ne even for less active stars (IRobrade et al. 
20081 ). As alluded to above, the vast majority of spectroscopically observed logLx < 29.1 
stars, including 7r 3 Ori, have Ne/O ~ 0.4. On ly a few truly inactive stars (the Sun and 
a Cen AB) have Ne/O w 0.2 dWood et alJboii ). 



There ar e a co uple characteristics of the low-FIP elements that are worthy of note. 
Maggio et al.l (120111 ) found that the coronal Ni abundance of the mid-F dwarf r Boo was 
anomalously high compared to other low-FIP elements, being a factor of 2 — 3 higher. We 
do not see this behavior for our mid-F dwarf, 7r 3 Ori. The Ni abundance may be a little high 
in Figure 5, but it is consistent with the other low-FIP elements within the error bars. As 
for Mg and Si, in the past we had noted a tendency for Mg a bundances to be highe r than Si 
abundances in the coronae of many moderately active stars (IWood fe Linskyl 120101 ) . This is 
not the ca se for 7r 3 Ori, t hough , which if anything shows the opposite behavior, consistent 
with what iMaggio et al.l (120111 ) found for r Boo. 
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Figure 5a shows the coronal abundance measurements relative to Fe. But using the 
log [Fe/H] = —4.32 measurement from the line-to-continuum analysis, we can compute ab- 
solute abundances as well, and in Figure 5b the absolute abundances are plotted, relative to 
stellar photospheric abundances as before. This figure casts doubt on whether the coronal 
abundance fractionation pattern of it 3 Ori is truly solar-like or not. 

Based on Figure 5a, the 7r 3 Ori coronal abundances were described above as very solar- 
like, because the Fi, ias = —0.55 measurement implies a factor of 4 underabundance of high- 
FIP elements relative to low-FIP elements in the corona, similar to the situation in the solar 
corona. But in the solar corona this is really due to an enhancement of low-FIP elements, as 
opposed to a depletion of high-FIP elements, meaning that for the Sun the high-FIP elements 
would be near 0.0 in Figure 5b and the low-FIP elements would be near 0.6. In contrast, 
Figure 5b would imply that for ir 3 Ori, the Ff, ias = —0.55 result is due to fractionation of 
both low-FIP and high-FIP elements, and is in fact mostly due to a depletion of high-FIP 
elements. 

T his is actually consiste nt with what we have found before for other moderately active 
stars (IWood fe Linskyl 120061 ) . The question is whether this means that the fractionation 
behavior of these stars really is significantly different from the Sun, or whether there is a 
systematic error in the line-to-continuum analysis that leads to systematic underestimates of 
log [Fe/H]. The history of solar coronal abundance measurements provides reason for caution 
with regards to the accuracy of the line-to-continuum analysis, as many such solar analyses 
have found the data more consistent with a high-F IP depletion than a low-FIP enhancement 
( TVeck fe Parkinson! Il98ll ; iFludra fe Schmelzl 119951 ). contradicting other evidence, including 
direct solar wind measurements, that suggest low -FIP enhancement is the primary effect 
( Ivon Steiger et al.lll995t iFeldman fe Lamina 120001 ). 



One concern is the possible presence of weak metal lines not listed in atomic line 
databases. These lines could collectively produce a significant pseudo-continuum in X-ray 
spectra, leading to overestimates of the strength of the real continuum, corresponding to 
an underestimate of log [Fe/H]. There is hope that modern databases such as CHIANTI 
with more extensive line lists have greatly reduced this problem, but it is difficult to quan- 
tify what the systematic uncertainties really are. Another potential issue is the effect of 
non-Maxwellian distributions on the line-to-continuum analysis. Non-Maxwellian distribu- 
tions will affect the ana lysis of both the line and continuum fluxes, but in different ways 
JPudik et al.lboill boij ). 



All these issues provide justification for relying primarily on relative, rather than ab- 
solute abundances. The latter introduces the uncertain systematic errors involved in the 
line-to-continuum analysis, while the former does not. Nevertheless, when using relative 



- 12 - 



abundances, and quantities derived from them such as F^as, it is important to keep in mind 
that the issue of absolute abundances has been left unresolved. A negative F bias value, for 
example, could indicate either a low-FIP enhancement or a high-FIP depletion, and two stars 
with identical F bias values could in principle be different in that respect. That being said, 
it is hard to imagine how the rather tight FBST relation shown in Figure 6 (and discussed 
in detail in the next section) could exist without some consistency regarding how F bias re- 
lates to absol ute abundances. F or example, the FBST relation suggests M4 V stars all have 
Fuas ~ 0.45 (jWood et al.l 120121 ) . It is hard to imagine why these stars would have different 
degrees of low-FIP depletions and/or high-FIP enhancements, but always calibrated so that 
Fbias « 0.45. 



3. REASSESSING THE CORONAL ABUNDANCE PROBLEM 



3.1. 7T 3 Ori and the FBST Relation 



As described in section 1, one of the main goals of the ir 3 Ori analysis is to see whether 
the FBST relation extends to F spectral types. T he answer is p rovid ed in Figure 6. A full 
list of the stars plotted in the figure is provided by IWood et al.l ( 120121 ). but we now add the 
7r 3 Ori data point. 

Based on the G, K, and M stars in the figure, tc 3 Ori would be expected to have 
Fbias ~ —0.8, potentially the strongest solar-like FIP effect ever observed. Instead it has 
Fbias — —0.55, comparable to the Sun and other early G stars. There are three possible 
interpretations for the higher-than-expected Fuas measurement for tt 3 Ori: 



1. Simply taken at face value, the 7r 3 Ori measurement indicates that the FBST relation 
flattens towards spectral types earlier than G. Observations of early F dwarfs would be 
very useful to verify this flattening. If real, the flattening of the FBST relation could 
be interpreted as an interesting manifestation of the narrowing convection zone. 

2. Although somewhat higher than expected, it could be argued that the tt 3 Ori mea- 
surement is still close enough to the expected F bias m —0.8 value to be considered 
consistent with a near-linear FBST relation, considering the amount of scatter that 
one sees in Figure 6. 



3. Perhaps ir 3 Ori is actually too active to precisely follow the FBST relation. It is known 
that stars with logLx > 29. 1 lie significantly above the FBST relation ( IWood et al. 
2012h . With logLx = 28.96 ( jSchmitt fc Liefke! l2004h . tt 3 Ori is close to this limit. If 
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the Lx limit is actually spectral type dependent, and lower for F dwarfs, ir 3 Ori would 
be expected to lie somewhat above the FBST curve. 



In section 1, there were two F dwarfs discussed that are known from previously published 
work to be inconsistent with the FBST relatio n: Pr ocyon (F5 IV- V) and r Boo (F7 V). Based 
on measurements from iRaassen et al.l (120021 ) and iMaggio et al.l (120111 ). we compute values 
of Fu a s = 0.12 and Fu as = —0.17 for Procyon and r Boo, respectively. These values are 
plotted in Figure 6, explicitly showing just how inconsistent these two stars are with both 
7r 3 Ori and the FBST relation. 

As for Procyon, the problem is that it is not a pure main sequence star, as its "IV- 
V" luminosity class and large 2.05 R radius demonstrate ( IKervella et al.l I2004J ) . Not only 
does Procyon have coronal abundances very different from tt 3 Ori, but it has a much cooler 
corona as well, with an EM peak near logT = 6.2 and no significant EM above logT = 6.4 
(IRaassen et al.ll2002l ). Cool coronae are associated with less active stars, but Procyon's X-ray 
luminosity is log Lx = 28.51, only a factor of 3 lower than n 3 Ori, and comparable to many of 
the stars in the Figure 6 sam ple with hotter coronae. For example, e Eri has a coronal X-ray 
luminosity of log Lx = 28.33 (jSchmitt fc Liefkdl2004l ). lower than Procyon, but it has an EM 
distribution similar to it 3 Ori (see Figure 2). We believe that Procyon's inconsistency with 
it 3 Ori emphasizes just how sensitive coronal properties can be to luminosity class. Procyon 
simply cannot be considered to be a main sequence star, and it is dramatically inconsistent 
with the FBST relation as a consequence. 

As for r Boo, the inconsistency of r Boo with 7r 3 Ori strengthens the argument made 
in section 1 that r Boo's high Fb ias value is due to one of its two companions. Either the 
M2 V companion, r Boo B, is contributing substantial flux to the XMM exposure, leading 
to erroneous abundance measurements, or the close-in giant planet, r Boo b, is somehow 
affecting the coronal abundances. The latter possibility is easily the more interesting of the 
two, as it would add a new twist to the current hot topic of star-planet interaction. 



3.2. What Determines Coronal Abundances? 

At this point it is worthwhile to summarize what we know about how coronal abundances 
depend on stellar properties. The stellar properties that we know affect Fb ias are: 

1. Spectral Type. This is the FBST relation illustrated in Figure 6. There are numerous 
basic stellar, photospheric, and convection zone properties that are known to be tightly 
correlated with spectral type for main sequence stars (stellar radius, temperature, 
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mass, surface gravity, convection zone depth, etc.). The existence of the FBST relation 
demonstrates that somehow one or more of these characteristics also determines coronal 
abundances. It is worth noting that a spectral type de pendence is also known to exist 



for T Tauri stars and very active main sequence stars (IGiidel et all 120071 ; I Wood et al. 



2012|). 



Activity and/or Rapid Rotation. The FBST relation seems to work very well for the 
vast majority of main sequence stars. However, there is clearly an activity dependence 
as well. This is most evident for extremely active st ars with IorLy > 29.1, which all 



lie significantly above the FBST relation in Figur e 6 (IWood et al.ll2012f) . This activity 



dependence is also clearly seen in Figure 37 of iGiidell (120041 ) . An indication that 



activity dependence of coronal abundances extends to less active stars is the difference 



in Ne/O seen 1 


)etween 


(Robrade et al. 


2008). 



For the activity extremes, it is possible that some aspect of 
high activity (e.g., high magnetic field strength) is affecting coronal abundances. But 
it is also possible that the underlying cause of the extreme activity, rapid rotation, is 
the true determining factor. Perhaps very rapid rotation affects photospheric and/or 
convection zone characteristics, thereby allowing coronal abundances to be altered in 
the same manner as spectral type affects abundances in the FBST relation. 

3. Luminosity Class. The best example of this is the ir 3 Ori/Procyon dichotomy discussed 
in section 3.1. There are evolved stars other than Procyon th at would lie aboye the 



FBST relation in Figure 6 , but most are part o f actiy e binaries (IBrinkman et all 12001 



Huenemoerder et all l200ll ; lAudard et all l200ll . 120031 ) . For such stars, it is not clear 



whether the discrepancy from the FBST relation is truly due to luminosity class, or 
whether it is instead an effect of high activity or in some cases particularly close 
binarity. 



Is there a common theme that can be found among these abundance determinants? One 
avenue to explore concerns photospheric oscillation properties, which are known to depend 
on both spectral type and luminosity class (#1 and #3 above), and may depend on rotation 
and activity as well (#2). If coronal abundance fractionation is cont rolled by Alfven and 



other MHD waves traveling through the chromosphere, as suggested by iLamingi (120041 . 12009 



2012h . the photospheric oscillation characteristics might be expected to be crucial, since 
these oscillations are likely sources of many of the chromospheric waves. 

Oscillation fre quencies increase and amplitudes decrease towards later type stars on 



the main sequence (IKjeldsen fc Bedding! 1 199 
F bias values 



20 111 ; IChaplin et al.ll2008l ). which have higher 
This would imply a potential connection between high Fbi as and oscillations 
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with high frequency and/or low amplitude. Unfortunately, this hypothetical connection is 
not consistent with the ir 3 Ori/Procyon dichotomy. Osci llation frequencies decrease and 



amplitudes increase as stars evolve off the main sequence ( IKjeldsen fc Bedding Il995l 12011 



Chaplin et al.ll2008l ). so the aforementioned main sequence oscillation/F 6ias connection would 
predict lower Fu as for Procyon compared with 7r 3 Ori. This is inconsistent with what we 
observe, so photospheric oscillation properties cannot be the sole determinant of coronal 
abundances. 



3.3. A Possible Theoretical Framework for Resolving the Coronal Abundance 

Problem 



Turning to theoretical considerations, we outline a framew ork for t r ying t o inte rpret 
the observed corona l abun dance variations in solar-like stars. lLamingi (120091 . 120121 ) and 
Rakowski fc Laming! (120121 ) argue that the Alfven waves giving rise to the ponderomotive 
force that generates the solar FIP effect most plausibly have a coronal origin. This is because 
waves in resonance with a coronal loop having period r = 2L/uva, where L is the coronal 
loop length and v a is the coronal Alfven speed, are required to produce the observed coronal 
depletion of He, and the only way to guarantee su ch a matching is if processes in the corona 
(e.g., nanoflares) excite normal modes of the loop. IWood et al.l (120121 ) present a "toy model" 
of the inverse FIP effect, where fast mode waves propagating up through the chromosphere 
from below are reflected back down again in the region of the chromosphere where the Alfven 
speed is increasing with height. This wave population gives a ponderomotive force directed 
downwards, which competes with the upwards ponderomotive force induced by the coronal 
Alfven waves propagating downwards and reflecting back up. We argue that at the F and G 
dwarf end of the FBST relation, the positive FIP effect from coronal Alfven waves dominates, 
while at the M dwarf end the inverse FIP effect due to initially upward propagating fast mode 
waves dominates, with a transition between the two somewhere in the K dwarfs. 

Why should this be so? We discuss first the coronal Alfven waves. These are assumed 
to ar i se from coronal nanoflare reconnec tion events. Several authors (e.g.. lLongcope fc Tarr 
2012 ; Kigure et al. 2010 ; Sturrock! Il999l ) have argued that at least some, and maybe even 
most, of the magnetic energy released in reconnection should be converted to kinetic energy 
or waves in the surrounding plasma. It is plausible that such wave generation explains why 
surveys to find localize d hot plasma as evi dence of nanoflare reconnection have generally 



been unsuccessful (e.g., I Warren et al.l 1201 lh . Instead, energy goes from magnetic field to 



waves, and is thus gradually dissipated as heat throughout the corona, and not quickly and 
locally as might have been expected. 
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Drake et al.l (120061 ) studied the efficiency of electron acceleration in reconnection with 
the ambient plasma (the electron plasma /5 e = 87rn e k-BT e /B 2 , to be more precise). They 
found maximum energy input to electrons at (3 e = 0, with reduced electron heating at higher 
(3 e , or lower magnetic field. We suggest that at the left hand side of the FBST relation, 
coronal reconnection primarily generates Aflven waves that end up causing positive FIP 
fractionation when they reflect from the chromosphere. As one moves to the right, to later 
spectral type, coronal reconnection puts more energy into electrons, and less into waves. 
Consequently, the positive FIP effect diminishes. 

Conversely, it is also possible that the amplitudes of fast mode waves penetrating the 
chromosphere from below and reflecting back down again may increase toward later spectral 
types, thus increasing t he inverse FIP effect. Even though p-mode amplit udes decrease 
towards later type stars ( IKjeldsen &: Beddindll995l . l2011t IChaplin et al.l 120081 ). it is possible 
that later type stars could be more efficient at converting the p-mode oscillations to fast 
mode MHD waves. Different regions of the Sun are diff erent in this respect. It is known that 
sunspots are sinks of p-mode energy (see lBraunlll995l . and references therein), most likely 
through mode conversion or resonant absorption to Alfven or fast mode waves. There are 
many cases of sunspots, active regions, and fla res observed to have reduced FIP effects (e.g. 
Feldman et al.lll990l ; iPhillips et al.lll994l 120061 ). although it is useful to note that an inverse 
FIP effect has never been reported in any solar observation. 

An important factor that could act to restrict inverse FIP effects to later spectral types 
is the expansion o f the magnetic field within the chromosphere. In their "toy model," 
Wood et al.l (120121 ) found that the density and magnetic field scale heights, Hp and Hb, 
were required to satisfy \Hd\ < \Hb\ /Q to yield an inverse FIP effect. Although more 
realistic models are required to refine this relation, we note that Hp should be smaller 
for later spectral types due to higher surface gravities, allowing the relation to be met more 
readily if H B does not also systematically decrease with spectral type. It is also worth noting 
that this relation will be met more readily for extremely active stars with magnetic fields 
so crowded that the fields cannot expand in the chromosphere, which effectively means an 
increase in Hb- This could potentially explain the tendency towards high Fb ias for extremely 
active stars (dependence #2 in section 3.2). 

The 7r 3 Ori/Procyon dichotomy (see section 3.1) remains difficult to explain within this 
framework. Procyon would have larger Hp than 7r 3 Ori, which would bias the star towards 
lower F bias , not higher (as observed), assuming magnetic field expansion is the central issue. 
If Hb increases by much more than Hp. as one moves from the main sequence to more evolved 
stars, that would reverse this behavior, but we know of no physical reason to expect this to 
be the case. 
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Perhaps Procyon has a fundamentally different coronal heating mechanism than the 
stars that follow the FBST relation. We have argued in section 2.4 that the commonly 
observed emission measure spike at logT ~ 6.6 would be inconsistent with a Type II spicule 
model of coronal heating and mass supply, if interpreted in terms of a thermal stability 
argument. But Procyon does not possess this peak (see section 3.1), so a spicule-based 
coronal heating picture might work for Procyon. Such a model could also be consistent with 
Procyon's lack of a FIP effect, because the spicule flow speed up through the chromosphere 
would conceivably be too fast to allow fractionation to occur. In this picture, the lower 
surface gravity of Procyon would allow upward propagating spicules to reach higher altitude 
and presumably higher temperature than on main sequence stars, potentially allowing the 
spicules to dominate the rather cool coronal emission observed from the star. Drake et al. 
(1995, their subsection 6.2) previously argued that the absence of FIP effect in the corona of 
Procyon is due to the effect of "unresolved fine structures" extending to higher temperatures 
there than they do on the Sun. 

These speculations provide a framework for trying to understand how coronal abun- 
dances depend on stellar parameters. We hope this discussion can serve to focus attention 
on the critical observational and theoretical points required to subject it to further scrutiny. 

We would like to thank Harry Warren for providing us with the solar active region EM 
distribution in Figure 2, and the referee Manuel Giidel for helpful comments. Support for 
this work was provided by NASA through ATP award NNH11AQ23I and Chandra Award 
Number GO1-12012Z issued by the Chandra X-ray Center (CXC). 
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Table 1. Chandra Line Measurements 



Ion A rest (A) logT Counts Flux (1(T 5 cm" 2 s" 1 ) 



Si XIII 


6.648 


6 


.99 


142.0 ± 30.3 


4.14 ±0.88 


Si XIII 


6.688 


6 


.99 






Si XIII 


6.740 


6 


.99 






Mg XII 


8.419 


7 


11 


56.9 ±20.2 


1.97 ±0.70 


Mg XII 


8.425 


7 


11 






Mg XI 


9.169 


6 


.80 


164.2 ± 27.3 


6.63 ± 1.10 


Mg XI 


9.231 


6 


.80 






Mg XI 


9.314 


6 


.79 






Nc X 


12.132 


6 


.87 


191.7 ±23.8 


8.89 ± 1.10 


Ne X 


12.138 


6 


.87 






Fc XVII 


12.264 


6 


.62 


117.3 ± 21.4 


5.47 ± 1.00 


Fe XXI 


12.285 


6 


.98 






Ni XIX 


12.435 


6 


.69 


74.2 ± 20.1 


3.48 ±0.94 


Ne IX 


13.447 


6 


.58 


314.3 ± 36.4 


14.41 ± 1.67 


Ne IX 


13.553 


6 


.58 






Ne IX 


13.699 


6 


.58 


98.2 ± 22.9 


4.51 ± 1.05 


Fc XVII 


13.823 


6 


.60 


173.2 ± 27.4 


7.97 ± 1.26 


Ni XIX 


14.043 


6 


.68 


39.1 ± 18.9 


1.80 ±0.87 


Ni XIX 


14.077 


6 


.68 






Fe XVIII 


14.203 


6 


.74 


292.1 ± 33.8 


13.47 ± 1.56 


Fc XVIII 


14.208 


6 


.74 






Fc XVII 


15.015 


6 


.59 


990.1 ±44.4 


44.83 ±2.01 


Fc XVII 


15.262 


6 


.59 


481.6 ± 34.7 


21.66 ± 1.56 


O VIII 


15.176 


6 


.65 






Fe XIX 


15.298 


6 


.83 






O VIII 


16.006 


6 


.63 


432.2 ± 34.2 


19.77 ± 1.56 


Fc XVIII 


16.005 


6 


.73 






Fe XVIII 


16.072 


6 


.73 






Fc XVII 


16.778 


6 


.58 


650.1 ± 36.4 


29.20 ± 1.63 


Fe XVII 


17.053 


6 


.58 


1428.3 ±47.8 


76.28 ±2.55 


Fe XVII 


17.098 


6 


.58 






O VII 


18.627 


6 


.34 


36.7 ± 17.0 


1.93 ±0.89 


O VIII 


18.967 


6 


.59 


624.5 ± 34.8 


32.58 ± 1.82 


O VIII 


18.973 


6 


.59 






O VII 


21.602 


6 


.32 


152.2 ± 24.5 


12.23 ± 1.97 


O VII 


21.807 


6 


.32 


36.5 ± 19.1 


2.98 ± 1.56 


O VII 


22.101 


6 


.31 


131.0 ± 23.9 


10.75 ± 1.96 


N VII 


24.779 


6 


.43 


47.5 ± 19.7 


3.97 ± 1.65 


N VII 


24.785 


6 


.43 






C VI 


33.734 


6 


.24 


114.0 ±21.2 


12.78 ±2.38 


C VI 


33.740 


6 


.24 






S XIII 


35.667 


6 


.43 


43.6 ± 18.6 


5.09 ±2.17 


Si XI 


43.763 


6 


.25 


94.4 ±20.9 


5.31 ± 1.18 


Si XII 


44.019 


6 


.44 


112.6 ± 21.4 


6.29 ± 1.20 


Si XII 


44.165 


6 


,44 


191.1 ± 23.8 


10.67 ± 1.33 


Si XII 


45.521 


6 


.44 


34.6 ± 19.1 


1.97 ± 1.09 
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Table 1 — Continued 



Ion A res t (A) logT Counts Flux (10~ 5 cm" 2 s" 1 ) 



Si XII 


45.691 


6 


.44 


qo 7 + 23 


U.iJU _l_ 1.01 


Fe XVI 


46.661 


5 


.43 


1 1 6 7 + 23 6 

J--LVJ. 1 _1_ • J . U 


6.77 ± 1.37 


Fe XVI 


46.718 


5 


.43 






Si XI 


49.222 


Q 


24 


<jj.u m j. y . o 


4 ^ -I- i ^ 
1.00 m ± . 00 


Fe XVI 


50.361 


5 


43 


1^1 1 4- 24 4 


1 ^ 81 + 2 ^ 

U.Ol _1_ £t>%J\J 


Fe XVI 




Q 


43 




q (?n _|_ 1 04 


Si X 


50.524 


Q 


.15 






Si X 


50.691 




15 




4 qq _i_ 9 91 


Si XI 


52.298 


5 


24 


14.9 ± 14.1 


2 1 ^ 4- 2 D4 


Fe XVI 


54.127 




43 


ftp. c _|_ on k 
uu.u m zu.u 


q nfi 4- 9 7Q 


Fe XVI 


54.710 




43 


i i 7 7 _|_ oo i 


1 q «7 _i_ 9 K7 


Mg X 


57.876 


Q 


22 


«4 Q _|_ 1 Q « 

u^t. o m j. y . u 


ft ne: _i_ 1 04 
u.uu m ±.01 


Mg X 


57.920 


5 


.22 






Fe XV 


59.405 


5 


32 


A ft _|_ 1 q 4 

uu . u m ± y . i 


d 4Q + 1 Q7 


Fe XVI 


62.872 




.42 


41 q _i_ 1 « 4 


c; 71 _|_ 9 97 


Mg X 


63.152 


6 


.21 


1 q + 1 3 f) 


9 04 _i_ 1 qq 


Mg X 


63.295 




21 


01 .u in ±u . ^ 


c; 70 _)_ 2 AO 
0. t u in ^ .*±y 


re A.V1 


Do. /II 


O. 




oo.y ± lo.y 


1 A 7*7 I fir 
ID. 1 I ± O.DO 


Fe XVI 


66.249 


6 


.42 


173.2 ± 26.1 


26.63 ±4.01 


Fe XVI 


66.357 


6 


.42 






Fe XV 


69.682 


6 


.32 


94.3 ±20.2 


9.97 ±2.14 


Fe XV 


69.941 


6 


.32 


22.3 ± 14.9 


2.38 ± 1.59 


Fe XV 


69.987 


6 


.32 






Fe XV 


70.054 


6 


.32 






Fe XV 


73.472 


6 


.32 


45.4 ± 17.9 


5.39 ±2.13 


Fe XVI 


76.497 


6 


.42 


42.2 ± 18.1 


5.52 ± 2.37 


Ne VIII 


88.082 


5 


.96 


50.2 ± 19.7 


8.46 ± 3.32 


Ne VIII 


88.120 


5. 


.96 






Fe XVIII 


93.923 


6 


.68 


303.0 ±36.8 


54.48 ±6.62 


Fe XIX 


101.550 


6 


.80 


29.8 ±20.1 


6.66 ± 4.49 


Fe XVIII 


103.937 


6 


.68 


82.0 ±29.1 


18.86 ±6.69 


Fe XIX 


108.355 


6 


.79 


133.1 ± 30.5 


30.85 ± 7.07 


Fe XX 


121.845 


6 


.88 


82.9 ±28.8 


23.14 ±8.04 


Fe XXI 


128.752 


6 


.95 


31.9 ± 19.9 


13.27 ±8.28 


Fe XX 


132.840 


6 


.88 


88.2 ± 32.0 


36.44 ± 13.22 


Fe XXIII 


132.906 


7. 


.12 






Fe XXII 


135.755 


7. 


.03 


38.9 ±24.1 


15.73 ±9.75 


Fe IX 


171.073 


5 


.95 


33.0 ±22.6 


46.92 ± 32.13 


Fe X 


174.534 


6 


.04 


33.9 ±22.5 


125.82 ±83.51 
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Table 2. Elemental Abundances 





Stellar 


Stellar 


Solar 


Abundance 


Corona 


Photosphere* 


Photosphere 13 


Fe/Fe 


1.5 


1.0 


1.0 


log [Fe/H] 


-4.32 


-4.50 


-4.50 


log[C/Fe] 


0.701°;^ 


1.12 


0.93 


log [N/Fe] 


-o-03±g : ig 


(0.49) 


0.33 


log [O/Fe] 


n 70+0.03 

u - ' °-0.04 


1.35 


1.19 


log [Ne/Fe] 


0-30^ 


(0.96) 


0.43 


log [Mg/Fe] 


0-081™ 


0.08 


0.10 


log [Si/Fe] 


n 16+ - 06 

u - iu -0.07 


0.07 


0.01 


log[S/Fe] 


-n i^+ - 16 

u - io -0.64 


(-0.22) 


-0.38 


log [Ni/Fe] 


i O1+0.12 
1 - zl -0.25 


-1.33 


-1.28 



a From Allende Prieto et al. (2004), except for values 
in parentheses, which are assumed rather than directly 
measured (see text). 

b From Asplund et al. (2009). 
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120 
Wavelength 

Fig. 1. — The Chandra LETGS spectrum of tc 3 Ori, rebirmed by a factor of 3 to improve 
S/N. For wavelengths above 35 A, the spectrum is also smoothed for the sake of appearance. 
The red line is a synthetic spectrum computed using the emission measure distribution in 
Figure 2, and the green lines indicate the contributions of higher spectral orders (2-5) to the 
model spectrum. 
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Fig. 2. — Coronal emission measure distribution derived for 7r 3 Ori, with 90% confidence 
error bars. This is compared with distributions for e Eri (dot-dashed line) and £ Boo A 
(dashed line), which were also measured from Chandra/ LETGS data (Wood & Linsky 2006, 
2010). Also shown for comparison is an average emission measure distribution for solar active 
regions (dotted line) studied by Warren et al. (2012). The solar distribution is normalized 
so that its peak matches that of 7r 3 Ori. 
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Fig. 3. — A synthetic line-plus-continuum spectrum (red line) showing the best fit to the 
highly smoothed Chandra 7r 3 Ori spectrum (gray line). The continuum contribution (blue 
line) to the synthetic spectrum assumes an absolute Fe abundance of log[Fe/H] = —4.32. 
The green line shows the contributions of higher spectral orders (2-5) to the total spectrum. 
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Fig. 4. — Illustration of four simulations used to assess the reliability of EM computations by 
PINTofALE. Dotted lines in each panels are emission measure distributions used to compute 
line fluxes, which are then fed into PINTofALE with the same relative uncertainties as 
the lines measured for 7r 3 Ori. Solid lines show the resulting EM distributions, with 90% 
confidence intervals. In all cases, these agree reasonably well with the actual distribution. 
For Simulation 1, the initial seed distribution is the ir 3 Ori distribution from Fig. 2. This 
distribution is also used for the other three simulations, but it is modified in the logT = 
6.5 — 6.7 range (marked by dashed lines) to explore PINTofALE's ability to reproduce EM 
peaks at these temperatures. 
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Fig. 5. — (a) The coronal abundance ratios of elements relative to Fe ([X/Fe]) are normalized 
to the photospheric abundance ratios ([X/Fe]*) and plotted versus first ionization potential 
(i.e., FIP), in eV. Error bars are 90% confidence intervals derived from the emission measure 
analysis. The vertical dotted line separates low-FIP elements from high-FIP elements. The 
horizontal dot-dashed line is the average relative abundance of the four principal high-FIP 
elements (C, N, O, Ne), which is the Fi, ias quantity used in Fig. 6. (b) Absolute coronal abun- 
dances of 7r 3 Ori relative to the photosphere, based on the log [Fe/H] = —4.32 measurement 
from the line-to-continuum ratio analysis (see Fig. 3). 
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Fig. 6. — A plot of FIP bias (Fb ias ) versus spectral type for the sample of main sequence 
stars from Wood et al. (2012), but with the addition of our new 7r 3 Ori measurement. 
Measurements are also shown for two additional F stars, Procyon and r Boo, which are 
clearly inconsistent with 7r 3 Ori and the general FlP-bias/spectral-type (FBST) relation. 



